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Star Formation Rates in Cooling Flow Clusters: A UV Pilot 
Study with Archival XMM-Newton Optical Monitor Data 

A.K. Hicks 1 and R. Mushotzky 2 

ABSTRACT 

We have analyzed XMM-Newton Optical Monitor (OM) UV (180-400 nm) 
data for a sample of 33 galaxies. 30 are cluster member galaxies, and nine of 
these are central cluster galaxies (CCGs) in cooling flow clusters having mass 
deposition rates which span a range of 8 — 525 M© yr _1 . By comparing the ratio 
of UV to 2MASS J band fluxes, we find a significant UV excess in many, but 
not all, cooling flow CCGs, a finding consistent with the outcome of previous 
studies based on optical imaging data (McNamara & O’Connell 1989; Cardiel, 
Gorgas, & Aragon-Salamanca 1998; Crawford et al. 1999). This UV excess is a 
direct indication of the presence of young massive stars, and therefore recent star 
formation, in these galaxies. Using the Starburst99 spectral energy distribution 
(SED) model of continuous star formation over a 900 Myr period, we derive star 
formation rates of 0.2 — 219 M© yr _1 for the cooling flow sample. For 2/3 of this 
sample it is possible to equate Chandra/XMM cooling flow mass deposition rates 
with UV inferred star formation rates, for a combination of starburst lifetime and 
IMF slope. This is a pilot study of the well populated XMM UV cluster archive 
and a more extensive follow up study is currently underway. 

Subject headings: cooling flows — galaxies:clusters:general — galaxiesrelliptical and 
lenticular, cD — galaxies: stellar content — stars: formation — ultraviolet : stars 

1. Introduction 

The ’mystery’ of cooling flows (Fabian et al. 1984; Sarazin 1986; Fabian 1994) has baffled 
astronomers for over 20 years. Though the X-ray temperature often drops by a factor of three 
within the central 100 kpc (Allen et al. 2001), there is seldom enough mass at temperatures 
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cooler than ~ 2 keV to be consistent with constant pressure cooling flow models. Recent 
XMM and Chandra results show that the distribution of temperature vs. emission measure 
also does not match the prediction of these models (Peterson et al. 2003; Voigt & Fabian 
2004). While the XMM and Chandra results have reduced the implied cooling rate in most 
clusters there is still an apparent discrepancy of an order or magnitude or more between the 
X-ray cooling rates and the observed star formation rates (SFR) in many clusters. Contrary 
to conventional wisdom that there is no star formation at all in the central cluster galaxies, 
a number of contemporary studies have found significant star formation (Melniek, Gopal- 
Krishna, & Terlevich 1997; Smith et al. 1997; Hansen, Jorgensen, & Norgaard-Nielsen 1998; 
Mittaz et al. 2001; McNamara et al. 2004a; O’Dea et al. 2004; McNamara 2004b). However 
we are still far from a complete reconciliation between X-ray mass deposition rates and 
alternatively derived star formation rates in the majority of cooling flow clusters. 

While many cluster cooling flows have been well observed in the radio and with emission 
lin e studies (Heckman et al. 1989, and others) there have been very few published studies 
of the UV properties of cooling flow clusters. In this study we use the readily available 
archival XMM Optical Monitor UV observations to derive the UV luminosities of a moderate 
sample of clusters. Since the dominant contribution to UV flux comes from short-lived main 
sequence stars, and the total number of these stars is directly proportional to the star 
formation rate (Kennicutt 1998) the UV is one of the prime routes to understanding star 
formation. We attempt to address and quantify the connection between cooling flows and 
star formation in our pilot sample. Unless otherwise noted, this paper assumes a cosmology 
of H 0 = 70 km s -1 Mpc -1 , Om = 0.3, and Qa == 0.7 


2. The XMM-Newton Optical Monitor 

The XMM Optical/UV Monitor Telescope is a modified 30cm Ritchey-Chretien tele- 
scope mounted alongside the X-ray mirror modules (Mason et al. 2001) which allows simul- 
taneous X-ray and UV/optical observations of XMM targets within the central 17 square 
arcminute region of the X-ray field of view. The 256 x 256 pixels have an angular resolution 
of ~ 1 " , The sensitivity limit of the Optical Monitor (OM) is 24 mag. This study primarily 
utilizes the UVW1 filter (220-400 nm), and the UVW2 filter (180-260 nm), and the UVM2 
filter (200-280 nm) for calibration purposes. Previous XMM OM studies of clusters have 
been limited to the detection of UV mission in Abell 1795 (Mittaz et al. 2001). 
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3. Sample and Observations 

Our pilot sample consists of 33 galaxies and eleven clusters. Three of them (M49, M60, 
and M86) were selected because they possess IUE photometric information (Burstein et al. 
1988) and were used to compare OM to IUE fluxes. 21 galaxies in the IC 1860 and Coma 
fields were used to calibrate the relationship between UV and IR luminosities in passively 
evolving early type and cD galaxies. Two UVW1 observations of Coma were included, one 
centered on the cD galaxies, and one which is off-center and contains numerous early-type 
galaxy detections within the viewing field. 

The remaining nine clusters in our sample were chosen because of significant X-ray de- 
termined mass deposition rates, general overlap with the samples of McNamara & O’Connell 
(1989) and Cardiel, Gorgas, & Aragon-Salamanca (1998), moderate to low redshift, inclu- 
sion in the 2MASS galaxy sample, and having archival XMM UVW1 data. The majority of 
clusters in our sample are low redshift objects (z < 0.1), but we have also included three 
moderate redshift clusters (z ~ 0.25) which are known to possess very strong cooling flows. 

Figure 1 shows one of the cooling flow clusters in our sample, with XMM-Newton UVW1 
contours overlaid. Table 1 lists the objects in our current sample, their redshifts, relevant 
filters, and OM exposure times. 


4. Surface Brightness 


We fitted the flux in concentric circles around the emission peak of each galaxy to a 
radial surface brightness profile of the form: 


/ „ 2 \ 1 / 2 — 3/3 

Hr) = Jo (l + ^) + 1 b (1) 

Where Jo is the normalization, tq is the core radius, and Ib is a constant representing the 
intensity contribution of the background. 

The surface brightness profiles are well fit by these models, cf. Figure 2. The fit 
parameters are in Table 2. 
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5. Photometry 

We retrieved pipeline-processed UV data from the XMM-Newton Science Archive and 
combined observations with the same pointing. Flux determinations for our calibration 
sample used a circular aperture of 7" radius, a good approximation of the IUE aperture 
(Burstein et al. 1988), and the size of the fixed 2MASS photometric aperture. Photon 
counts in the aperture were corrected for background using a 35" radius background region 
taken from an outlying area of the chip in the same observation. Final fluxes were determined 
by employing the OM counts-to-flux conversions, and errors were assessed via -\/N photon 
statistics. 

The results from UVW2 and UVM2 observations of M49, M60, and M86 were then 
compared with the IUE determined flux value of each object (Burstein et al. 1988) for the 
bandpass of the appropriate filter. Our flux measurements are consistent with theirs to within 
< 35%, errors which are in keeping with those of Vassiliadis et al (1996) and Koratkar et al. 
(1997). Two of the three galaxies have lower OM fluxes than IUE fluxes, indicating that the 
Optical Monitor does not have a red leak. 

We then determined UVW1 fluxes for the remainder of our sample (see Table 1), for 
a 7" radius region, in the method described above. The UVW1 data were corrected for 
average galactic extinctions (Cardelli, Clayton, & Mathis 1989). J band extinction values 
were retrieved from the NASA/IPAC Extragalactic Database (NED). 


6. Calibration 

The low redshift (z ~ 0.023) calibration sample has 23 galaxies in IC 1860 and Coma 
with 2MASS J band photometric measurements. Two of these galaxies exhibited high UV 
fluxes with respect to the rest of the calibration sample, and were described as star forming 
in a FOCA study Donas, Millard, & Laget (1995) and thus are not appropriate for inclusion 
in our UV/IR luminosity calibration of passively evolving early type galaxies, and were 
removed from the remainder of this study. 

A least squares fit between the UVW1 and J band luminosities for the remaining 21 
early type galaxies is well described by 


Luv = C\ x Lj^ 2 


( 2 ) 


Where Luv is the UVW1 luminosity and Li R is the 2MASS J band luminosity, inside 



-5- 


7" radius apertures with C\ — 0.090 ± 0.007 and C' 2 = 0.708 ± 0.007 respectively (Figure 
3). We use this correlation to predict the UV light from non-starforming, passive galaxies 
assuming that the UV light from such objects is produced by old evolved stars coming from 
the same population that produces the IR light. 


7. Photometric Results 

The UV luminosity excesses of the cD galaxies in our nine cooling flow clusters were 
calculated by subtracting the expected UV light due to old stars (based on the fit obtained 
in Section 6) from the measured UV luminosity of each cluster. Figure 4 shows the J band 
luminosity of each galaxy in our sample plotted against its UV/IR luminosity ratio. We have 
not attempted to estimate and correct for internal dust absorption, and thus these excesses 
provide a lower limit on the UV emission in these clusters. 

Because we are analyzing data from a single filter, it is difficult to break the degeneracy 
between age and mass of the young stellar populations. We folded through the UVWl filter 
Starburst99 (Leitherer 1999) redshifted models corresponding to continuous star formation 
over a a variety of timescales with solar abundance (Z = 0.020) and IMF powerlaw indices 
of both 2.35 and 3.3 to establish the relation between OM UV counting rates and star 
formation rate. The measured UV luminosity excesses of the clusters were then used to 
estimate per-year star formation rates (Table 3) that were then compared to X-ray derived 
mass deposition rates from Voigt & Fabian (2004), Peterson et al. (2003), and R. Mushotzky 
(2005, private communication). Also in Table 3 are estimates of the duration of continuous 
star formation, and best fitting IMF powerlaw index. These values were derived by plotting 
the relationship between star formation duration and rate for each given UV excess and IMF 
slope, then finding the set of parameters (slope and duration) at which X-ray determined 
mass deposition rates intercepted those curves. The 900 Myr models should be closest to 
the average lifetime of a cooling flow (Allen et al. 2001). 


8. Summary and Discussion 

We find a significant UV luminosity excess in the central galaxies of most cooling flow 
clusters. This implies that star formation has recently been or is currently occurring most 
likely as a direct result of radiative cooling of the intracluster medium (ICM). These con- 
clusions are qualitatively consistent with the findings of previous optical imaging studies of 
star formation in cooling flow clusters (McNamara & O’Connell 1989; Cardiel, Gorgas, & 
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Aragon- Salamanca 1998; Crawford et al. 1999). 

In the absence of multiwavelength spectral information, it is difficult to constrain a 
particular model of the star formation history of these galaxies. However we note that for 
2/3 of our sample it is possible to find some combination of the parameters of IMF slope and 
starburst lifetime which allow us to equate X-ray determined mass deposition rates with UV 
inferred star formation rates. 

In none of the OM, Chandra or HST images is there evidence of a significant excess 
central flux, therefore we expect the AGN UV contribution to our signal is small. In addition, 
we have investigated the amount of UV flux which might come from the cooling flow itself, 
and have found it to be much smaller than measured UV excesses. Since the calibration 
sample consists entirely of low redshift early type galaxies, there is a possible systematic 
error in our calculations of the SFRs of the three moderate redshift (z ~ 0.25) clusters in 
our sample, due to evolution in the UV luminosity of the old stellar population (Brown et 
al. 2003) and/or uncertainties in the redshift evolution of the star forming component. In 
future work we will attempt to calibrate a higher redshift passive E galaxy sample. 

While there is significant UV flux in 2/3 of our cooling flow sample, 1/3 of the objects 
with significant cooling rates from modern data do not have detectable UV excesses, and 
the ratio of UV light to X-ray cooling rate (Figure 5) shows a wide range. The origins of 
this scatter are unclear and a larger UV sample, with direct comparison of X-ray and radio 
imaging and IR data, is necessary to get a clearer picture of the origin of the UV light. The 
discovery of copious UV emission from several cooling flow clusters confirms the previous 
work and shows that a significant fraction of the cooling gas in clusters can indeed form 
stars. 

Our star formation estimates are very basic, being based solely on a simple conversion 
from UV light to SFR (Kennicutt 1998). We have bypassed issues of the age of the star 
formation, the effects of reddening, the slope of the IMF and metallicity effects. The fact 
that we can, roughly, obtain similar values of X-ray inferred cooling rate and UV inferred 
star formation rate indicates that, given the above degeneracies, more and better data (such 
as could be obtained by Galex or HST spectroscopy) is necessary to determine precise star 
formation rates. 

This paper is a pilot study and we are currently undertaking an expansion of this 
project toward higher redshift, utilizing more of the OM UV filters and including the many 
additional clusters that remain in the archive. Through this expansion we hope to achieve a 
more thorough characterization of recent star formation in cooling flow clusters of galaxies. 
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(220-400 
exposure 

image, and have count values of 14, 25, 44.75, and 80. The large splotch of UV contours to 
the lower right of the cD galaxy is due to scattered light. 
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Fig. 2. — Surface Brightness Profile. Radial surface brightness profile of Abell 1795 
made from UVW1 (220-400 nm) data accumulated in ~ 1" annular bins. A solid line traces 
the best (3 model fit. The horizontal dotted line represents the best fitting background value. 
The best fitting (3 model has a core radius of 4.09 " and beta value of 0.703 
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J Band Luminosity vs. UV Luminosity 



Fig. 3. — UV vs. IR luminosity for non-cooling flow cluster galaxies in our sample. The line 
indicates the best fitting relationship between UV and IR luminosity for passively evolving 
and non-cooling flow cD galaxies, represented by squares and triangles, respectively. Dia- 
monds indicate the two star forming galaxies (Donas, Millard, & Laget 1995) which were 
excluded from the fit. A relationship between UV and IR luminosity is expected here, be- 
cause passively evolving cluster galaxies have similar spectral energy distributions (SEDs). 
The similar redshifts of IC1860 and Coma allow us to fit this relationship without making 
distance corrections to the data. 
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Fig. 4.— UV/IR luminosity ratio of selected galaxies vs. their 2MASS J band luminosities. 
All luminosities were calculated in 1" radius regions, and corrected for galactic extinction. 
Diamonds indicate cooling flow cluster cDs, all of which are at low redshift (z < 0.1) except 
for the three points farthest to the right, which have redshifts of ~ 0.25. The line designates 
the best fitting relationship between UV and IR luminosity for passively evolving and non- 
cooling flow cD galaxies (Section 6). The high UV/IR ratio of the cooling flow clusters is a 
strong indication of the existence of hot young stars and thus a direct indication of copious 
recent star formation. 
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Excess UV luminosity vs. M 



X-ray Mass Deposition Rate [M 0 /yr] 


Fig. 5. — Excess UV luminosity vs. M. Luminosity excesses were determined by subtracting 
the best fit ratio of UV light due to old stars (Figure 2) from the measured UV luminosity 
of each cluster. X-ray mass deposition rates were taken from recent papers (Peterson et al. 
2003; Voigt & Fabian 2004). 
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Table 1. XMM-Newton Optical Monitor Observations 


Object 

z 

Filter 

Exposure [s] 

IC 1860 

0.0229 

UVW1 

5000 

PKS 0745-191 

0.1028 

UVW1 

3999 

Hydra A 

0.0538 

UVW1 

3000 

Zw3146 

0.2906 

UVW1 

4000 

NGC 4406 (M86) 

-0.0008 

UVW1 

3600 



IJVM2 

3500 

NGC 4472 (M49) 

0.0033 

UVW1 

2000 



UVM2 

1000 

NGC 4649 (M60) 

0.0037 

UVW1 

2499 



UVW2 

2900 

Coma (center) 

0.0231 

UVW1 

3900 

Coma (outskirts) 


UVW1 

1498 

Abell 1795 

0.0633 

UVW1 

2499 

Abell 1835 

0.2520 

UVW1 

2000 

MS 1455.0+2232 

0.2578 

UVW1 

2000 

Abell 2029 

0.0779 

UVW1 

1200 

Abell 2052 

0.0345 

UVW1 

2000 

MKW3s 

0.0453 

UVW1 

2002 
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Table 2. Beta-Model Fits 


Object 

tfcore ["] 

P 

/ o a 

/ B a 

Reduced x 2 

PKS 0745-191 

2.83 ±0 ' 04 

0 766±°- 011 

625. 3 ±12,4 

481.0 ±L6 

0.691 

Hydra A 

3.449 ±a001 

0.8387 ±0 - 0002 

4829. 1 ±2 - 4 

356. 59 ±006 

4.200 

Zw3146 

6.21 ±001 

2.054 ±0 - 007 

1119. 2 ±4-3 

413.2 ±0 - 4 

0.947 

Abell 1795 

4 09±°- 01 

0.703 ±o -° 02 

2334. 8 ±13-4 

336. 6 ±0 - 3 

0.946 

Abell 1835 

3.65 ±0 ' 41 

1.253 ±ai65 

2880.7 ±153 - 0 

349 6±i- 8 

0.514 

MS1455. 0+2232 

4 99±°.°3 

1.006 ±o -° 07 

446 9± 3 -5 

337.8 ±0 - 4 

0.657 

Abell 2029 

3.68 ±0 - 02 

0.499 ±0 - 001 

613 2 ± 3 - 5 

449. 7 ±a6 

1.101 

Abell 2052 

1 .68 ±a28 

0 402±°- 022 

1593. 1 ±149 ' 6 

426. 2 ±7 - 4 

1.912 

MKW3s 

1 ,28 ±0 - 33 

O.352 ±0 - 022 

1018.0 ±138 ' 9 

398. 6 ±71 

1.110 


a Surface brightness I in units of 10 17 ergs sec 1 cm 2 arcsec 2 
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Table 3. 


Cluster 

z 

UV Excess 
[10 43 erg/s] 

^2.35 a 

[Mq/yA 

M3.3 a 
[M G /yr] 

M (X-ray) 
Wq/ yx] 

7 sf 

[Myr] 

a 

PKS 0745-191 

0.1028 

11 94±°.7 

129 ±7 

237 ±l3 

317 + || b 

200 

3.3 

Hydra A 

0.05384 

0.88 ±0 -° 2 

9.5 ±0 - 2 

17.5 ±0 - 4 

3+6 b 

900+ 

2.35 

Zw3146 

0.2906 

8.47 ±0 - 4 

91± 4 

168* 8 

525±lg 

11 

3.3 

Abell 1795 

0.06326 

0.85 ±a ° 2 

9.1 ±0 - 2 

16.7 ±0 - 4 

30^ b 

100 

3.3 

Abell 1835 

0.251982 

11.41 ±cu 

123 ±5 

226 ±9 

34±|| b 

900+ 

2.35 

MS 1455.0+2232 

0.2578 

3.36 ±0 - 4 

36 ±4 

66 ±8 

200 + 7q 

150 

3.3 

Abell 2029 

0.07795 

0.1 8 ±0 ’ 05 

y±0.5 

3 2 ±°.9 

3oiS b 

7 

3.3 

Abell 2052 

0.03446 

0.1 9 ±0 * 007 

0.19 ±o -° 8 

0.4 ±ai 

ioo± 2 ° c 

... 


MKW3s 

0.04531 

o.oi ±001 

0.1 ±ai 

0.2 ±a2 

45«° c 




a Assuming constant star formation over 900 Myr. The subscript indicates the IMF power- 
law slope of the model used. 

b Voigt & Fabian (2004). 

c Peterson et al. (2003). 



